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We analyze the event-by-event fluctuations of mean transverse-momentum measured recently by the PHENIX 
and STAR Collaborations at RHIC. We argue that the observed scaling of strength of dynamical fluctuations with 
the inverse number of particles can be naturally explained by formation of multiparticle clusters. 
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I. INTRODUCTION 

This paper is based on our recent analysis fjl of the event- 
by-event fluctuations of the mean transverse-momentum mea- 
sured by PHENIX H and STAR J3 0] Collaborations at 
RHIC. The discussed new data has extended our knowl- 
edge of the phenomena studied intensely in recent years 
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the central issues in this field is whether the pj fluctuations at 
intermediate momenta result from jets 12211 - In this paper 
we argue that the observed scaling of strength of dynamical 
fluctuations may be naturally explained by formation of mul- 
tiparticle clusters. From our point of view, the explanation of 
the data in terms of quenched jets represents only one possible 
option, other physical mechanisms leading to clustering also 
explain the data. 

Our reasoning is based on the following three observations: 
the mean and the variance of the inclusive momentum dis- 
tribution are practically constant at low centrality parameters 
(c = 0-30%), 
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it) the variance of the average momenta for the mixed events 
is practically equal to the variance of the inclusive distribution 
divided by the average multiplicity, 
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the extracted value of Od yn is large at the expected scale pro- 
vided by the variance of pj, which indicates that the clusters 
should contain at least several particles in order to combina- 
torically enhance the magnitude to the observed level. This 
point is supported by the numerical calculations of the values 
of (Jjj coming from the resonance decays and from thermal 
clusters in statistical models of heavy-ion collisions, which 
show that these values are very small | ll 12311 - 



n. SIMPLE WAY OF DOING STATISTICAL ANALYSIS 

In this Section we present elementary statistical consider- 
ations which form the basis of our treatment of the fluctua- 
tions. We use the notation where the letter p denotes \pr\, 
Pi is the value of p for the zth particle, and M = YJ!=\ Pi/n 
is the mean transverse momentum in an event of multiplic- 



TABLE I: Analysis of the event-by-event fluctuations in the trans- 
verse momentum. Upper rows: the PHENIX experimental data at 
^Jsnn = 130 GeV |2]; middle rows: the mixed-event results; bottom 
rows: our way of looking at the data. One observes that to a good 
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approximation a M 
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: o],/{n) and a% n = (a 2 , - a^ mlx ) ~ l/(n>. 
Except for the first two rows, all values are given in MeV. The errors 
in the last row reflect the unknown round-off errors in the data of the 
upper and middle parts. 



finally, Hi) the difference of the experimental and mixed-event 
variances of mean pj, denoted as C7^ yn , scales as inverse mul- 
tiplicity 
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°dyn = °M ~ °M 
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As mentioned above, an explanation of this scaling can be 
provided by clustering in the expansion velocity. Moreover, 
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ity n. Consider events of multiplicity n and transverse mo- 
menta pi,p2, ■ ■ ■ ,p„. The multiplicity n and the momenta 
are varying randomly from event to event. The probability 
density of occurrence of a given momentum configuration is 
P(n)p„(p\ ,...,/?„), where P(n) is the multiplicity distribution 
and p n (pi, • • • ,Pn) is the conditional probability distribution 
of occurrence of pi,...,p„ in accepted events, provided we 
have the multiplicity n. Note that in general p depends func- 
tionally on n, as indicated by the subscript. The normalization 



£p(n) = l, dpi...dp„p„(pi,...,p„) = l. (4) 

n 

The marginal probability densities are defined as 

Pn~ k \pu---,P„-k) = J dp n - k+ l...dp n p n (pi,...,p n ), (5) 

with k = 1 ,...,«— 1 . These are also normalized to unity, as 
follows from Eq. (@}. Since the number of arguments distin- 

(n—k) 

guishes the marginal distributions p„ , in the following we 
drop the superscript (n — k). Further, we introduce the follow- 
ing definitions 



dpPn(p)p, 



(6) 
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(P)n = 

var n (p) = ]dp 9n (p)(p-tp) 

COV n (pup 2 ) = J dpidp 2 (pi - (p)n){P2- (p)„)Pn(pUP2)- 

The subscript n indicates that the averaging is taken in samples 
of multiplicity n. We note in passing that the commonly used 
inclusive distributions are related to the marginal probability 
distributions in the following way: 



PinM =£f(") / dp\...dp„Y,§(x- Pi)Pn(Pl,---,Pn) 
n •> i=\ 

= £«P(«)p„W, (7) 

n 

9m{x,y) 

r n 

/ dp\...dp n £ 5{x-pi)5(y-pj)p„(pu---,Pn) 

n 7 i,;= ijy; 

= "£n(n-l)P(n)p n (x,y), (8) 

n 

which are normalized to (n) and (n(n — 1)), respectively. For 
the variable M = XH=i Pi/ n we f" nc l immediately 
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(M 2 



Y, p ( n ) / dpi...dp n Mp„(p u ...,p n ) 
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Y,P{n) / dpi...dp„M 2 p n (pi 
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FIG. 1: The cluster picture of correlations. Particles are grouped 
in iV c i clusters, containing on the average r particles. The particles 
within a cluster move at very similar collective velocities, indicated 
by arrows. 



III. ANALYSIS OF THE PHENIX DATA 

At first we discuss the PHENIX measurement |2] of the 
event-by-event fluctuations of the mean transverse momen- 
tum at y/SNN =130 GeV. We use the experimental fact that 
the variance of the momentum distribution and its mean are 
independent of centrality (in the considered range), which al- 
lows us to replace the quantities (p}„ by (M) and (p 2 )„ — {p) 2 
by Op at the average multiplicity, denoted as a 2 /„\. In this 



way we get 



J M 



P(n) 



P(n) 



£ COV„(pi,Pj) 
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In the mixed events, by construction, particles are not corre- 
lated, hence the covariance term in Eq. il Q vanishes and 



2,mix _ 2 ' 
°M - a p,(n) t 



Pin) 



, (12) 



where in the last equality we have used the fact that the distri- 
bution P(n) is narrow and expanded l/n = 1 /[(«) + (« — («))] 
to second order in (n — (n)). Comparison made in TableU see 
the 8th and 9th rows, shows that formula (I12> works at the 
1-2% level. In addition, since O p j„) is not altered by the event 
mixing procedure, subtracting dl 2I > from il It yields 



°dyn = L 



P(n) 



£ COV n (pi,Pj) 
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TT2 L COV {n) (pi, P j). 



(13) 
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IV. CLUSTER EXPLANATION 

The scaling of the dynamical fluctuations with («), see the 
last row of Tableland Eq. (0, may be regarded as a constraint 
restricting the behavior of the covariance term. For instance, 
one can immediately conclude that if all particles were corre- 
lated to each other, L" ;=1 /^cov^ B )(p,-,/?y) would be propor- 
tional to the number of pairs, and Od yn would not depend on 
(n) at large multiplicities. A natural explanation of the scaling 
(0 comes from the cluster model, depicted in Fig.^ The sys- 
tem is assumed to have N Q \ clusters, each containing (on the 
average) r particles. Below we keep r = const, for simplicity 
(a more general case is discussed in Ref. | The particles 
are correlated if and only if they belong to the same cluster, 
where the covariance per pair is 2cov*. The number of corre- 
lated pairs within a cluster is r(r— l)/2. Some particles may 
be unclustered, hence the ratio of clustered to all particles is 
(N c })r/(n) = a. If all particles are clustered then a = 1. With 
these assumptions Eq. (1131 becomes 



J dyn 



a(r- 1) 



(14) 



which complies to the scaling <|3j . An immediate conclusion 
here is that the ratio a cannot depend on («) (in the considered 
centrality range, c = — 30%) in order for the scaling to hold. 

The question now is whether we can use the above results 
to draw conclusions on effects of jets (minijets), which have 
been proposed as a possible explanation of the experimental 
data even at the considered soft momenta 12111 . Jets, when 
fragmenting, lead to clusters in the momentum space. The 

resulting full covariance from jets is then N c ij e tj(j — 1)t^> 
where /V c ij et is the number of clusters originating from jets, 
j is the average number of particles in the cluster, and 2 cov J 



is the average covariance per pair. The total number of par- 
ticles produced from jets is N c \j et j. On the other hand, the 
commonly accepted estimate of the dependence of /VcijetV' on 
centrality is accounted for by the nuclear modification factor 
Raa multiplied by the number of binary nucleon-nucleon col- 
lisons iVbin- Since Raa depends on the ratio («) / (n) pp , where 
(n) P p is the multiplicity in the proton-proton collisions, in a 
given pr bin one finds 



Ncljelj ~ RAANbin = 



N bin ~( n ), (15) 



Nbin(n/ P p 

which complies to the scaling of Eq. dl4> . 

V. CONCLUSIONS 

The observed scaling of the dynamical fluctuations of mean 
transverse-momentum may be naturally explained by the clus- 
ter model. We emphasize that the scaling follows from the 
presence of clusters only, and is insensitive to the nature of 
their physical origin as long as one imposes the condition 
N d ~{n). 

The explanation of the observed data in terms of quenched 
jets agrees with the cluster picture. However, the explanation 
of the centrality dependence of the pr fluctuations in terms of 
jets based solely on Eq. Jl 5i is not decisive: any mechanism 
leading to clusters would do, for instance the thermal clusters 
discussed in Ref. (j! or mechanisms of Ref. i24ll . Micro- 
scopic realistic estimates of the magnitude of cov J and j are 
necessary in that regard, including the interplay of the jets and 
medium. In short, the nature of clusters remains an open issue. 
For the current status of this program the reader is referred to 
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